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Abstract 

An effective Hamiltonian which could model the interaction between a tunnehng proton and 
the conduction electrons of a metal is investigated. A remarkably simple correlation between the 
motion of the TLS-atoui and an angular-momentum change of scattering electron is deduced, at 
the first-order Born level, by using a momentum-space representation with plane waves for initial 
and final states. It is shown that the angular average of the scattering amplitude-change at the 
Fermi surface depends solely on the difference of the first two phase shifts, for small-distance 
displacements of the heavy particle. For such a limit of displacement, and within a distorted- wave 
Born approximation for initial and final states, the change in the scattering amplitude is expressed 
via trigonometric functions of scattering phase shifts at the Fermi energy. The numerical value 
of this change is analyzed in the framework of a self-consistent screening description for impurity- 
embedding in a paramagnetic electron gas. In order to discuss the so-called antiabatic limit on 
the same footing, a comparison with matrix elements obtained by the potential-gradient of an 
unscreened Coulomb field is given as well. The coupling of the tunneling proton to a free-electron- 
like electron gas is in the typical range obtained, by ultrasound experiments for different metallic 
glasses, from scattering rates for a Korringa-type relaxation process. That coupling is too weak to 
be in the range required for realization of the two-channel Kondo effect. 



PACS numbers: 72.10.Fk, 72.15.Qni 



I. INTRODUCTION 



In the last decades hydrogens in metals deserved very extensive experimental studies and 
vast theoretical considerations. In crystalline solid the hydrogen (H) sits in a well-defined 
interatomic position like, e.g., in Pd or Pt. That is not the case in amorphous systems and 
at dislocations and other nonperiodic distortions. If the H has some more room between 
the host atoms its position may be not well-defined and it moves between two positions. 
Such systems are known as two-level systems (TLS) and they have been very extensively 
studied^. The coupling to the conduction electrons can result in an extra contribution to 
the electrical resistivity^. If the atom has two metastable positions the electron scattering 
amplitude in different angular momentum channels depends on the atomic position. The 
difference between these amplitudes is described in the literature by a coupling V^. That 
coupling contributes to the resistivity in a conventional way. 

There are, however, other couplings where electrons induce transitions between the two 
levels. Thus an assisted transition can be realized by the tunneling of the atom between 
the two positions^i^i^. The importance of that coupling, denoted by and V^, is highly 
debated^i^. These models consider different atoms with sizeable differences in their masses. 
The original suggestion is limited to small tunneling rate, while more intensive tunneling 
induces an essential split between the energies of the atomic eigenstates which reduces their 
roles. The model has attracted considerable interest as it was suggested that at low tem- 
perature exhibits non Fermi liquid behavior known as the two-channel Kondo (2CK) effect, 
which has its own theoretical interest^i^. 

Experimental supports come much more from the studies of point contacts than from 
direct measurement of the electrical resistivity^. The Cornell group made detailed suggestion 
how the observed zero-bias anomalies could be due to that 2CK effect which has been also 
highly debated^. Since that time the original model was modified^ by taking into account 
the actual electronic structure at the TLS. The possibility was also considered^ where the 
atom moves between the two positions via the next higher energy level of the atomic motion. 
These suggestions were aimed to increase the electron assisted amplitudes {V^, V"^) to make 
the 2CK more feasible. The enhancement of these couplings by the renormalization due to 
the conduction electrons is drived by V^. The estimated couplings were on the borderline, 
therefore, those should be studied in more details. It is crucial whether those coupling 
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strength could reach some certain regions which are very sensitive on the strength of V^. 

Recently it was suggested that if is a possible candidate^. Indeed, zero-bias anomalies 
have been observed in hydrogenated Pd point contact with considerable size^i. That system 
could deserve more extensive experimental studies. Such zero-bias anomalies were also 
observed in hydrogenated Pt point contact^. The idea has beeen also raised that other 
systems may contain some water with metastable positions of the hydrogen atoms. 

These recent developments justifies further studies of the coupling of a heavy tunneling 
particle to conduction electrons. The satisfactory realiable estimation of these coupling 
could be performed by considering certain atomic configuration with detailed knowledge 
of the electronic density of states in the region of the TLS. In other words, a material- 
specific description of the host could be based upon more realistic but much more complex 
specifications. Thus, a calculation of the coupling would be too ambitious in case where the 
structural surrounding would require more knowledge. 

The present paper is devoted to a minimal program. An embedded proton is considered 
in a double-well potential which determines the two possible positions of the heavy parti- 
cle. These positions are considered as stable ones. The main task is to take into account 
the screening, due to the Coulomb interaction and charge-response, in a satisfactory way. 
In this work the ground-state screening in the surrounding electron gas is described in a 
self-consistent Hartree-like manner, by using Kohn-Sham independent single-particle states 
whose occupation is prescribed by the Pauli principle. The change in the scattering ampli- 
tude, due to a potential-displacement, is expressed in our renormalized perturbation theory 
via matrix-elements taken between stationary scattering eigenstates of the external Kohn- 
Sham field, instead of commonly applied matrix-elements between unperturbed plane-wave 
states of a system without impurity. Concretely, in the present work we shall apply the 
distorted wave Born approximation^^ for scattering characteristics. 

The paper is organized as follows. The next section. Sec. II, is devoted to the details 
of our theory and a discussion of the plane-wave-based, first-order Born approximation for 
matrix-elements is also given. These results, obtained by using screened potentials, are 
compared with those based on an unscreened Coulomb field. The last section. Sec. Ill, 
contains a short summary and an outlook for further possible developments. The Appendix 
summarizes the commonly applied basic elements of previous motivating works. 
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FIG. 1: Illustration of the two-level system. The shadowed circle represents the screened proton 
in the central position and the two others, separated by d, the TLS. The dashed lines correspond 
to the spherical wave centered around the origin. The inset shows the potential of the TLS. 

II. THEORY AND RESULTS 

As we motivated in the Introduction, an atom is considered which has two stable positions 
in metallic matrix due to a double potential well. The positions are symmetric to a central 
point and they are at z = respectively. This is the standard picture to a theoretically 
active field with considerable experimental relevances; see Figure 1 for illustration. 

The main goal of the present contribution is to determine the change in the electron 
scattering amplitude in different angular-momentum channels when the atom is moved out 
from the central point to one of the two positions, but the set of the electron wave functions 
is still centered at the origin. These positions, at 2; = ±| as the illustrative Figure shows, 
are described by a pseudospin = ±1. 

The atom is embedded in the degenerate electron gas of the metallic target. The gen- 
eral form of the Hamiltonian is H = Hq + Hi, where the diagonal Hq matrix stands for 
independent electrons in stationary eigenstates of a self-consistent external field V{r) 



in which the are energy eigenvalues of bound and scattering eigenstates. The a' and a. 




(1) 
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create and annihilate these eigenstates of Eq.(l) of spin a. 

In order to show the difference of our distorted- wave Born method from the commonly 
applied plane-wave {pw) approximation, and thus provide a clear phenomenology to under- 
standing, we start by assuming a weak potential-energy and write h[^^^ in this unperturbed- 
state representation for initial and final states as 

= ^'^T.^pMT.^Un^^^ (2) 

q p 

with the corresponding operators for creation and annihilation of plane-wave states. Here 
we have q = k2 — ki, and thus Vp^{q) is given by 

= I dve-'^-^^V{v). (3) 

Since the perturbation AV(r), due to the shift in the atomic position, is real we have 
the [^^(q)]* = q) character. In other words, the h[^^^ is Hermitian. The change 

in the scattering amplitude is AFpw{q) = [m/(27r/i^)]\^^(q). The extended continuous 
system can be recovered simply by using the prescription (l/V^)Xlq ~^ I c^q/(27r)'^. The 
q-summation must be performed at the Fermi level to respect the Pauli's principle, and 
thus the determinant-character— of the state vector of Hq. This will result in, by using 
the q = 2kFS\.n{d/2) exact expression for the momentum change in elastic scattering, an 
averaging (see below) over the scattering angle. From this point of view, our procedure 
is similar to the one^*^ used in the theory of conventional superconductivity to define a 
repulsive (electron-electron) pseudopotential, from a spherical V{q) at the Fermi surface. 

One can calculate the change [AF(k2, ki)] in the scattering amplitude beyond the usual 
Born approximation, by using expansions for properly defined^ initial and final states in 
terms of normalized spherical harmonics. Furthermore, it was suggested^ that one can define 
a Hermitian ifi, beyond the weak-coupling limit, in the following way 

^1 = ^ 5Z ^k.k. «k,., (4) 

kik2 

where V^kaki — //3(^2)Ki/3(^i' ^2)/a(^i) IS a suitable^ coupling, mediated now, in our 
case, by a Al^(r) perturbation. The intruduced^ fim{k) = i'' y'^Yimik) functions are spheri- 
cal harmonics, and a and (3 run over a properly chosen set (see below) of angular momentum 
indices {l,m). The present calculation is carried out in several steps. 
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(i) Calculations of the screened potential around the proton and the phase shifts. The 
charge is sitting in the central position where its screened field V{r) and a complete set 
of one-electron wave functions of occupied states are determined in a self-consistent way. 
This is achieved, in practice, by applying the Kohn-Sham method of density-functional 
theory (DFT), which reduces the complicated many-body problem of the inhomogeneous 
electron gas (in the presence of a charge) to a single particle problero^^. The calculations are 
performed for a grand-canonical system, i.e., at a fixed chemical potential, applying local- 
density approximation (LDA) for the exchange-correlation potential. The single-particle 
potential energy V{r) has a simple form in this approximation 

V{r) = -^+ [ dh' + Au^^Mr)], (5) 

r J \r — r I 

in which An(r) is the screening density. The many-body term Az/^-cfn] is expressed via an 
input exchange-correlation chemical potential (/i^c) as Au^dn] = fj.xc{no + An) — f^xci^o), in 
order to have a vanishing effective potential energy at infinity. 

For a given density (no) of the screening environment, and depending on the magnitude 
of the attractive embedded charge, the total density consists of bound and scattering eigen- 
states. The n(r) = no + An(r) total density, the basic variable of DFT, is constructed out 
by summing over doubly occupied bound and scattering-like states 

^±(r) = AKA;)/2K^,^)>i;,(^)//.n(^), (6) 

Im 

in which Ri{k,r) are self-consistent solutions of the radial Kohn-Sham equations with V{r) 
at {hky/2m scattering energy, and Ai{k) = e^*'^'*^'^-* where 5i(k) is the phase shift. The 
above continuum states are normalized on the fc-scale, and thus the scattering (sc) part of 
the induced density comes from an integral over the Fermi-Dirac distribution function 

Auscir) = - 5^ (2/ + 1) / dkP [RUk,r)-jf{kr)] . (7) 

^ 1=0 

The total-screening condition, An drr'^An{r) = Z, implies the Friedel sum rule of scat- 
tering phase shifts in one-electron mean-field treatments 

oo 2 °° 

Z = - V (2/ + 1) 5i{kF) - - V (2/ + 1) 5i{0) + m. (8) 

TT ^ — ' 71 ^ — ' 

1=0 1=0 
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FIG. 2: Self-consistently determined radial wave functions, Ro{kF,r) and Ri{kF,r), as a function 
of the radial distance r in atomic units (a.u.). Solid and dashed curves refer, respectively, to the 
I = and / = 1 components. The Wigner-Seitz parameter of the host system is = 2.5 (a.u.). 

The rule is satisfied, of course, at numerical self-consistency of iterations. Here A^;, denotes 
the number of occupied bound states. In a ground-state calculation, on which the present 
work is based, the last two terms cancel each other according to Levinson's theorem^. At 
metallic densities, already the first few phase-shifts provide a very accurate approximation. 
For example, at the = 2.5 value of the density parameter one has 5o(^f) = 1-2213 and 
^i(^f) = 0.0894 in the case of a proton. This effective can characterize the mobile part 
of the electron fluid of a Pd target, and satisfies^ a necessary condition appropriate for 
correctly describing ground state properties of defects. 

Illustrative results are exhibited in Figure 2 for the leading, I = and / = 1, ra- 
dial wave functions calculated at the = 2.5 value of the Wigner-Seitz parameter. The 
/ = component shows a Coulomb-like enhancement, ^y2n/kp at the origin r = 0, over 
the plane-wave-based unity. The deviations from the plane-wave components joikpr) and 
ji{kpr) are notable, as Figure 3 shows via the corresponding products of components. The 
RQ{kp,r)Ri{kF,r) product has its maximum at r = 0.7 with value 0.393, and it becomes 
zero at r=2.53. For higher r values it oscillates with decreasing amplitudes around zero due 
to (shifted) Friedel oscillations. As Figure 3 shows, the perturbative jo{kFr)ji{kFr) product 
is much more extended. Its maximum (0.26) is at about r=1.7, it becomes zero at about 
r=3.9, and has Friedel oscillations beyond this. Clearly, the attractive proton is a strong 
perturbation at metallic densities. 
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FIG. 3: Products of the leading partial waves as a function of r. The solid and dashed curves corre- 
spond, respectively, to the self-consistent Ro{kp,r)Ri{kF,r) and the perturbative Jo(^F^)ji(^F'")- 
The Wigner-Seitz parameter of the host system is = 2.5 (a.u.). 

As a digression we enumerate at this point few important facts which support the reha- 
bility of the above phase shifts and the radial Ri{kp, r) functions even in a real lattice where 
the structural (atomic) surrounding could have influence. First, even unit charges {Z = ±1) 
represent strong local perturbations. This is verified by experiments for slowing down of 
low-speed {v < vp) protons and antiprotons moving through paramagnetic metallic targets 
characterized by Svr^no = {pp/K}^, where pp = mvp = hkp. In the theoretical^^ energy loss 
per unit path length [{dE/dx) = mvvpnQatr{pp)] appears the transport-cross section 

4 °° 

(TtriPp) = (^p^jf^Y + ^) siv?[Si{pp) - 8l+i{pp)]. (9) 

The theoretical results, based on DFT phase shifts for protons^i and antiprotons^^, are in 
impressive agreement with {dE / dxYdata. obtained^^ at dedicated facilities of CERN and 
show, via the atr{Z = l,pp)/atr{Z = —l,pp) ratio (/?), a pronounced {R > 2) charge-sign 
effect which rules out^ the applicability of a simple Born-approximation. For completeness, 
for an antiproton one has 5o(^f)=— 0.7729 and Si{kp)=— 0.2003 at = 2.5. The illustra- 
tive Figures for the radial wave functions, and the above-outlined phase shifts values for 
embedded unit charges {Z = ±1) suggest [see the discussion at Eq.(24) also] the use of the 
so-called sp-approximation in order to get a reasonable estimation for our coupling constant. 
This will be, therefore, the practical approximation in the TLS problem. 
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Quite remarkably, a dimensionless parameter Koveriap{d), which expressea^^*^ the overlap 
of two many-body ground states with the same local potential at two sites separated by a 
short distance d, depends on the phase shifts similarly 

1 f 2d\'^ °° 

Koverlapid) = 3 ( ) 5Z + ^) sin^ipp) - (10) 
^ ^ ^ 1=0 

where = 27rh/pF is the de Broglie wavelength of an electron moving with vp velocity; 
(2d/XB) = {dhp/Ti). Related to positive muon (/i"*^) quantum diffusion in metallic targets, 
realistic phase shifts were already applied^^ to estimate this dimensionless parameter. 

Second, while the band-structure paradigm emphasizes the importance of lattice struc- 
ture, calculation^^ for the Pd — H system using a molecular-cluster model shows that the 
electronic properties of the impurity are dictated mainly by its local environment. The 
such-calculated on-top screening density, An(r = 0), is only slightly smaller than the sim- 
ple jellium-based result which tends^ to l/vr from above as no decreases. Furthermore, the 
molecular-cluster and pseudo-jellium calculations for Pd—H are in close agreement with each 
other for if-displacements up to 0.5ao from the equilibrium configuration. A self-consistent 
pseudopotential calculation^^ shows that the proton is screened on one atomic distance, and 
the hydrogen always has more charge around it than the Pd upto the Wigner-Seitz radius. 
Note, that the amount of electron-localization in iJ-screening could characterize, as was 
pointed out^ recently, the site-preference of hydrogen in new storage materials. 

We finish our supporting enumeration by associating the above /oca/-environment picture 
with the proposal of Hopfield^ for short-range properties when there is a change in potential 
due to moving an atom by a small distance which results in matrix-elements needed, in his 
case, to an estimation of electron-phonon coupling in transition-metal superconductivity. 
Namely, it was shown that when an angular momentum decomposition of electron wave 
functions is used, the matrix elements contain chiefly scatterings which change the angular 
momentum of the electron. Precisely, it is this character which is central in the context of 
scattering of electrons from two-level systems. The parity change of the angular momentum 
state (without altering the spin indices) of conduction electrons gives them an internal degree 
of freedom coupled to that of the impurity^i^. This internal degree is the background to 
establish an analogy with the usual (spin-related, magnetic) Kondo effect. 
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{a) The change of the potential for d ^ 0. The screened charge is moved to one of the 
positions cr^ = ±1 carrying the potential, which is taken rigid as the screening action is very 
fast^ compared to the infrared processes essential in some other problems. The electron 
eigenfunctions are also carried with the tunneling atom, but they are decomposed is terms 
of those have already been determined for the central position; see in point (i). Thus the 
new potentials are 

Vir Ta'^z) = V{r) ± a' Ar(r), (11) 
where the important (perturbative at small d) term is given by 

AT^(r) = V{r-U) - y(r). (12) 

A Taylor expansion in d results in a simple dipolar form in the leading order 

, , d dV(r) d yV(r) 

AV(r) = - — ^ = - cos^— 13) 
^ ^ 2 dr^ 2 dr ^ ^ 

We illustrate the behaviours of the self-consistent potential and its gradient in Figure 4. 
Concretely, the —rV{r) and r'^[dV{r)/dr] products are plotted by solid and dashed curves, 
respectively. The inset is devoted to the finer details of these important functions. The 
calculated results refer to = 2.5 for the density parameter. The remarkable Coulombic 
character of r'^[dV{r)/dr] at short distance is due to a compensating effect between the 
electrostatic (Hartree-term) screening and the local exchange-correlation term. 

In the following perturbation theories are applied in terms of AV{r). The matrix elements 
are calculated between the states with wave functions determined in the previous 

section (i) but with the plane wave functions are also presented. Namely, we use the distorted 
wave Born and the usual Born approximations, respectively. 

{iii) Matrix elements of the shifted potential between the original wave functions. In order 
to get a convenient, dimensionless^ coupling to characterize the effect of the perturbation 
on precalculated continuous states we shall use energy-normalization^i for these: 



The original matrix element at kp is multiplied, in such a way, by the density of states 
Po{Ef) = {kprn)/ {27i'^h'^) = 0.75no/EF per unit volume for a given spin evaluated at the 
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FIG. 4: Characteristics of the self-consistent potential obtained at = 2.5 for the screening of an 
embedded proton. The -rV{r) and r'^[dV{r)/dr] are plotted as a function of r, by solid and dashed 
curves, respectively. The inset is devoted to finer details for r > 3. Atomic units are used. 

Fermi energy. Notice, oncemore, that the AF{ki, change in the scattering amphtude 
[i.e., the matrix element of the perturbing AV^(r) between iP^{t) initial and final states^] 
should involve a m/{2Txfi^) prefactor. With Eq.(13) the angle-integration over VLf gives, by 
applying standard recurrence relation for cos 6 Yim{0,(t)), the following simple result 

= j d^f coseYi^{f)Y*,^,{f) = 

Only the /' = /+! and m' = m values are allowed, due to the dipolar character. 

Born approximation. It is instructive first to study the case of a weak, Vps{r) pseudo- 
potential in Eq.(13). Thus, we perform the radial integration in first-order Born approxi- 
mation, i.e., we apply plane- wave components for the free radial wave functions as 

h{l) = ^ drr'^^J^Ukr)j,^,{kr). (16) 

In this perturbative case the integration by parts and use of the following expression based 
on recurrence relations for Bessel functions 

^ [r'^ji{kir)ji+i{k2r)] = r'^^jK^ir) - fciJ?+i(A;2r)], (17) 

results in a remaining integral now with Vps{r). In such a way we can apply the standard 
definition^ of the first-order Born {B) phase-shift 



m 



(15) 
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^Fik) = drr'V,s{r)jKkr), (18) 

to obtain the following informative [see also Eq.(26), below] expression 

= [Sf{k)-5Uk)]. (19) 

We stress that this equation is valid, physically, only for small values of the phase shifts, i.e., 
when the distortion of the electron wave function by the central potential field is negligable. 

If a weak potential is given via its Fourier representation, as in a dielectric screening of 
the external field Vps{q) = Vextiq) / ^{q) , the Born phase shifts (we use atomic units here) are 



^F(k) = -^l dqqV,MF(^-lJ + l;l;^y (20) 

where F is the standard hypergeometric function. In our simple, so-called sp-approximation, 
one can get from Eq.(20) easily [see Eq.(24), also] the leading difference-term 

6^{k) - 5f{k) = -J^s I dqq' [-qVM]- (21) 

This g-space representation could be useful if the inverse Fourier-Hankel transfomation re- 
sults in a complicated function for Vps{r), and a fast estimation is needed. For a com- 
monly applied Yukawa-type potential, V^g{r) = —{Z/r)exp{—\r), the Born phase shifts 
from Eq.(18) are given by Legendre-functions [(5«(x)] of the second kind as 

Srik,\) = ^Qi{l + {\/kf/2). (22) 
It is important to note, in the present context, that a restricted Friedel sum 
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Z = - ^ (2/ + 1) 5r{k,, A) = - 4 V^M = 0), (23) 

1=0 

still holds despite the perturbative approximation z/A^ = Akp/TT, i.e., the screening param- 
eter corresponds to the quasiclassical Thomas-Fermi value. Physically, the normalization 
of the screening charge (calculated from the model Yukawa form via Poisson's equation) is 
satisfied, but its real-space distribution is not necessarily realistic. 

The first-order Born phase shift does not contain the multiple scattering effect in the 
central potential field. Formal applications of Eqs.(18)-(19) with strong potentials might 
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result in an uncontrollable estimation for the effect in question. However, the proposed 
correlation^ between the motion of the TLS-atom and an angular-momentum change is 
now transparent even at the first-order Born level, as Eq.(19) clearly shows. Note that 
this transparency of correlation is obtained (easily) when one applies^ spherical-harmonics- 
based expansions for the initial and final (plane waves in first-order Born approximation) 
states separately, i.e., by implementing Eq.(4). 

When, still in first-order Born approximation with a Vps{r) in Eq.(13), we implement 
Eqs.(2-3) with the standard spherical-harmonics expansion for e~^'^'^, i.e., without the men- 
tioned separation for initial and final states, one gets {Vpg = Vp^) for small enough d 



In this representation of the (perturbative) matrix-element a parity-change is not transpar- 
ent. But a partial integration and the use of the j'i{x) = jo{x) — {2/x)ji{x) recurrence relation 
gives —qVps{q) for the integral. The Fermi-surface average, (q^qdq) =^ [Poi^) — Pi{x)]dx 
with X = cos-i?, becomes, as Eq.(21) shows, a function of [5(f (/c^) — Si^kp)] solely. This 
conclusion agrees with Ref.[6]; only the / = and / = 1 harmonics are relevant. 

Beyond the first-order Born approximation. In a quite recent theoretical work-, which 
also rests on matrix element calculation with a potential-gradient between normalized s and 
p bound-states, the bare [Vc'(r) = —Z/r] Coulomb potential was applied to characterize TLS 
in a metallic matrix. Here, by using the 77 = Ze^m/ {kfi?) Sommerfeld parameter, we add 
the corresponding exact phase-shift difference for the case of a Coulomb field 



By using the accurate Ri{k, r) and Ri^i{k, r) radial wave functions to integration in Eq.(16) 
above [and not Rf{k,r) and Rf_^-^{k,r) to Eq.(18)] with gradient of the true V{r) behind 
these functions one obtains, based on earlier result^*^, the exact 

Ik{l) = sm[Si{k) - Si+i{k)] = [tan(5,(A;) - tan(5;+i(A;)] cos(5;(A;) cos5;+i(A;), (26) 

closed expression in terms of scattering phase shifts [see: Eqs.(9) and (10), also]. The exact 
result, which contains now the multiple scattering effect in the central potential field to 
all order, is a bounded function in contrast to Eq.(19). We illustrate, in Figure 5, the 




(24) 



5i{k) - %i(fc) = arctan 



V 



(25) 



(^ + 1) 



13 



0.4 








2 



3 



r(a.u.) 



FIG. 5: Illustration of the integrand of Eq.(16) obtained by using self-consistent result for the 
potential gradient and the two leading radial wave functions. The Wigner-Seitz parameter of the 
host system is = 2.5. Atomic units are used. 

argumentum-function r'^RQ{kp, r)Ri{kF, r)[dV{r)/dr] of Eq.(16) by using the self- consistent 
solutions at = 2.5. Fortunately, as noted above, we have already a simple, analytic result 
in Eq.(26) for the integral without further numerics. 

We add at this important point that with Eq.(25) and Eq.(26) one can get the 



expression in terms of the parameters Z, k, and /, for a Coulomb field; 1^(1) < 1. The 
plane-wave-based perturbative Coulomb limit, obtained with V{r) = —Ze^jr to Eq.(16) is 
surprisingly similar to Eq.(27). The difference is only, but this is crucial since /c^ ~ 1 in a.u., 
that there is not if term [see Eq.(21) also] in the denominator. Formal use of a finite Imax to 
Eq.(27) could mimic a screening-regularization and, together with Eq.(15), could thus allow 
fast analytical estimations. For the Coulomb case in the rj ^ oo limit (formally ^ ^ 0) one 
has [SQ^k) — Si{k)] [7r/2 — hve/Ze^], where Ve = hk/m is the electron velocity. Notice, 
parenthetically, that this difference is one of the pedagogical examples in physics which 
shows transparently the role of Planck's constant, "in action". 

After the above detailed analysis on scattering, we return now to the matrix element 
needed to a suitable^ coupling introduced in Eq.(4). First we give, by using the previous 
notations, a dimensionless expression related to the scattering-amplitude change 




V 



(27) 



[(/ + l)2 + r/2]V2' 
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A/^^(fcp) = (A;^/7r) AF„^^(A;^) = — I{l,m) h^l) A+^ikp) Aikp), (28) 

in which a and (3 refer to the {l,m) and (/ + l,m) values, respectively. Next, following 
earlier works^i^, we write a dimensionless form to the Hermitian Hi in Eq.(4) as 

^Po{EF)V:^{kF). (29) 

The obvious connection in the weak-perturbation limit is, based on the peculiarity of the 
first-order Born (B) approximation with plane waves, simply Vaif\^F) ~ ^f^p\^F)- 

This transparent connection between physical quantities at the Fermi surface suggests us 
to use, beyond the above weak-coupling limit, the vl^pi^kp) = Re ^f^pikp) extension. With 
our choice for boundary conditions to select initial and final states (involved in matrix- 
element calculation based on the distorted wave Born method) this seems to be the only 
logical step which preserves the important Hermitian character of Hi and reproduces the 
weak-coupling limit. In addition, the proposed extension is in harmony with standard text- 
book statement^ on the characterization of an energy shift of a particle interacting with a 
potential. Considering the above-deduced, remarkably simple rule on the true effect of an 
angle-averaging at the Fermi surface [see, at Eq.(24)], we write 

i?eA/o^oio(M = -}-^sin[6o{kF)-6i{kF)] cos[5o(M + ^i(M], (30) 

ImAf^^i n^ikp) = sin[5o(A;F) - Siikp)] sin[5o(A;i.) + (5i(/c^)] (31) 

for the important real (-Re) part, and the imaginary (Jm) part which is at least second-order 
in a weak-perturbation {Z — > 0) limit. 

Our physically-motivated extension to fix a value to the coupling of a TLS with con- 
duction electrons, remains within the framework of common knowledge: All stationary 
properties of metals which can be described in terms of scattering of conduction electrons 
based on an adiabatic picture, are periodic functions of the phase shift, except, curiously, 
the Friedel sum rule. We stress that in our model an electron merely sees a different scat- 
tering potential for each state of the configuration but has no effective internal-spin degree 
of freedom. The ijj^ states, for initial and final states, are needed^ to a perturbation theory 
in the continuum of a given potential field. In experiments for transport characteristics, we 
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have current-carrying electrons. Thus a standing-wave representation (which would refer to 
a different boundary condition described by physics, for example by a cavity) for the states 
involved in our matrix-element calculation is not reasonable. 

In the so-called unitary limit, So{kp) ~ 7r/2, where the effect of the self-consistent cen- 
tral V{r) field is strong [see. Figs. 2 and 3], the influence of the potential-shift [A\/(r)] 
becomes very small, i.e., "yooiol^^) small. Remarkably, this statement is in accord with 
the conclusion of Gyorffy on a renormalized limit (obtained via mapping to the partition 
function of a fictitious, auxiliary logarithmic gas) with contact interaction^. This shows 
that the ground-state charge-distribution around an embedded proton is, in fact, intact to 
small (~ d) extra fields. In a one-phase- shift approximation the maximal coupling would 
be at the 5o(^f) = vr/4 mathematical value. At this value the real and imaginary parts 
of the scattering amplitude are equal in magnitude. We shall return, briefly, to the simple 
contact-potential approximation in the Appendix. 

The expression in Eq.(28), which gives a strict linearity in d, is based on a leading- 
term expansion for the perturbation AV"(r) as Eq.(13) shows. In order to get a more 
detailed ^-dependence of coupling due to the displacement of a screened proton in the z- 
direction, we performed numerical volume-integrations with Eq.(13) and the dominating 
Ro^kp, r)YoQ{f)Ri{kF, r)YiQ{f) product. By introducing the simplifying notation 

K^{d) ^|A/oVo(M|Ai., (32) 

motivated by the need (see below) of explicit estimation for a physical observable, and 
performing the ^-integration via a variable-change, we have 

2\/3 rr+d/2 

K'{d) = ^ drRo{kF,r)Ri{kF,r) duuV{u)[r^ + {d/2Y - u% (33) 

" Jo J\r-d/2\ 

The numerical result for K^{d) is presented in Figure 6 by a solid curve. The dashed 
curve refers to the asymptotic expansion, which is linear in d. We note that the asymp- 
totic expansion provides a quite acceptable representation upto about d ^ 0.5, and even 
at c? = 1 the deviation from the numerical results is only about 25%. The somewhat sur- 
prising linearity can be explained, partly, by the fact that the second term (~ cos^ 6) of a 
formal Taylor series for the perturbation would give vanishing contribution, i.e., there is not 
quadratic, ci^-proportional, term at small d in the sp approximation. Beyond the physically 
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FIG. 6: Numerical result (solid curve) for K^[d), defined in Eq.(33), as a function of {d/2). The 
dashed curve refers to the asymptotic Taylor expansion which gives a linear dependence on d. The 
Wigner-Seitz parameter of the host system is = 2.5. Atomic units are used. 

reasonable d = 1 value, the K^(d) function grows more and more gradually and it has a 
maximum at about c? ^ 1.8. With an acceptable mathematical accuracy one may use the 
K^{d) ^ 0.52d (1 + 0.9d) exp{—0.9d) fit for < 3 in our sp-dominated problem. 

Finally, the magnitude of |A/qq ^0(^^7)1 can be calculated by using our phase shifts at 
density parameter = 2.5 [in this case kp = (97r/4)i/Vr, ~ 0.77], and assuming (see Fig. 
6) a size of the TLS d = 1 and then 



That value is moderate enough to ignore the multiple scattering due to AV (see Appendix B), 
and is in accordance^i^ with the estimated typical values for different metallic glassy systems 
in the intermediate coupling regions using ultrasound measurements. In our modelling of 
the amplitude-change due to a TLS, the absolute maximum at d ^ 1.8 does not provide 
more than an about 25% increase. A naive mixed approximation, in which the gradient of 
a bare Coulomb potential is weighted by our self-consistent Ro{kp,r)Ri{kF,r) product in 
Eq.(16), gives a sma//er numerical value than Eq.(26). This is due to, mathematically, an 
over-weighting of the negative part of the product in Fig. 3 by an unscreened-field gradient 
in the volume-integral. A careful discussion of the underlying adiabatic and anti-adiabatic 
pictures was given by Kagart^, by stressing the relevance of an adiabatic character due to 
the fast screening- act ion, in TLS-motion. 




(34) 
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III. SUMMARY AND OUTLOOK 



A dimensionless coupling constant that characterizes the effect of a potential-gradient 
perturbation on scattering eigenstates of a self-consistently treated embedded impurity is 
deduced for small values of the impurity displacement d in metallic electron gases. The result 
is expressed via bounded, trigonometric functions of scattering phase shift differences at the 
Fermi energy, quite similarly to a well-known overlap parameter Koveriap{d)- Beyond the 
leading-term expansion for the perturbation, the numerical results show that the coupling 
parameter has an almost linear (i-dependence upto the physically realistic d = Iqq value for 
the displacement of a screened proton in an electron gas with = 2.5. 

As we mentioned earlier, at the enumeration of supporting facts related to the applied self- 
consistent-field method, it is the parity change of the angular momentum state of conduction 
electrons which gives them an internal degree of freedom in the TLS problem. We derived 
this parity change in the present paper by using a distorted wave Born approximation, i.e., 
using those scattering wave functions to matrix-element calculation which are determined by 
the central field of a fixed, spherically screened proton. Beyond the applied distorted wave 
Born approximation, with an axially symmetric scattering potential, the exact description 
would lead to coupled radial equations^ in a partial wave expansion of the scattered wave at 
a given k. Particularly, the ci-dependence of the coupling beyond the present renormalized 
perturbation theory with precalculated spherical states, is an exciting problem. Application 
of internally consistent (axially symmetric) nonperturbative wave functions to the TLS 
problem needs a future study. The numerics could be based on a Z^FT-method developed 
recently— for axially symmetric potentials in an electron gas. 

It is generally true that a clean theory, like the present one based on scattering aspects, is 
an analysis of the properties of an idealized, hypothetical model. What we have tried to do is 
to give a self-contained mathematical treatment of a physically motivated model, hydrogen 
in an electron gas, that can demonstrate the strong local effect in an atomic displacement in a 
metallic target. A more detailed, i.e., material-specific description of the host could be based 
upon more realistic but much more complex physical specifications. For example, simple 
application of a specified, as a function of the hydrogen concentration, density of states at 
the corresponding Fermi energy of a real Pd — H system could change the numerical value 
of the dimensionless coupling, found here with a free-electron form pq{Ep) = 0.75no/ Ep. 
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Our self-consistent Kohn-Sham approximation incorporates the electron-electron interac- 
tion in the Hartree-like mean-field potential which acts for independent electrons. Beyond 
this approximation, and especially with bound states around the bare impurity, the treat- 
ment of electron-electron interaction (statistical and dynamical) has to be somewhat more 
sophisticated, since part of it is already included in the screening. One must avoid double 
counting, and consider self-interaction corrections. A strong interaction between the elec- 
trons, involved in the screening action, may well lead to the occurrence of localized magnetic 
moments associated with an embedded impurity atom. In such a cases the interplay between 
the orbital and spin degrees of freedom may become an important question. 

As we fixed in the present paper, in our model an electron merely sees a different scattering 
potential for each state of the (impurity) configuration but has no effective internal-spin 
degree of freedom. A nontrivial extension of the applied method, i.e., the distorted wave 
Born approximation, could be a problem where in addition to our displacement AV{r) 
there is an extra spin-orbit (so) coupling already at the level of our Hq via a scattering 
interaction: V{r) =^ V{r) + Vso{r), where Vso{r) ~ [r' -^V {r)L-(T]. As was demonstrated^ 
by Ballentine, the change in the scattering amplitude has a peculiar character in this case. 
Namely, there is no term with / = in the amplitude-change, because the operator L yields 
zero in that case. How a combined Al^(r) + Vsoi^) perturbation could change the statement 
obtained at Eq. (24) in a TLS problem might deserve a detailed study. 

The theoretical modeling of important observables based on transport-related experi- 
ments needs additional care. For example, a recent calculation shows the possibly important 
role of different charge-states in current-driven electromigration and backflow,— where the 
long-range characters of scattered waves are more important than in our present problem for 
a TLS coupling. We have discussed and emphasized the importance of physical boundary 
conditions (which appear in the fundamental Lippmann-Schwinger equation) for continu- 
ous states which are based on a second-order differential (Schrodinger) equation. If our 
displacing-atom were embedded into a system described by externally generated^ standing 
waves, one could use^ the e^^^^'^'^ ^1/ cos5/(/c) normalization-change and, thus, get for the 
right-hand-side of Eq.(26) only the difference of tangents. In this case, the coupling could 
enhance almost resonantly. Indeed, the mentioned STS study^ heralds, via spectroscopic 
informations, a common nonmagnetic effect with different impurities in standing-wave pat- 
terns. We speculate that the dissolution of hydrogen into Prf-electrodes can make a local 
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confinement-like effect for states involved in the conduct ance^i of Pd — H nanoj unctions. 

Finally, a way to consider electronic inhomogeneities can be an additional local-density 
approximation governed by the strong short-range distortion on which our model is based. 
Between Pd atoms of a real lattice there are ranges, measured from a lattice-atom position, 
where the density of states has enhanced local values. If the allowed tunneling occurs in 
such geometrical ranges, a mathematical averaging of our t;^[rs(r)] function over a certain 
range of rs{r) seems to be reasonable. The role of an almost ferromagnetic nature of a real 
Pd target (expressed via a Stoner-enhancement^ in the spin response function) might also 
deserve future considerations. These combined considerations could give further quantitative 
information to a field of considerable experimental relevance. 
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IV. APPENDIX A: CONNECTION WITH EARLIER THEORIES 

This Section is devoted to a comparison with previous woiksr-'^^ on modeling the coupling 
in Eq.(2), based on plane- wave states. In order to provide a clear phenomenology, we stress 
the point that the present theoretical description, which is based on calculation of matrix- 
elements of a dipolar potential-perturbation, also uses prefixed basis sets. Namely, plane- 
waves in the conventional and precalculated self-consistent ones in the distorted-wave Born 
approximation for continuous states. 

In the main text, at the details of the first-order Born approximation, we outlined the q- 
representation of the coupling. Thus the desired link, for example to Ref. [46] with a Yukawa 
potential, is easily obtained by using the scattering value of g = 2kpsm{'d/2) and the 

VM = = ^ E(2/ + l)5f (A^i.) Piicos^), (35) 

Born-representation to Eq.(24) for the V^^(q) quantity; P/(cos'/9) ~ Yi^(ki)Yim{k2) ■ 
Motivated by certain scattering-length arguments, the ai notations were adopted in the 
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mentioned earlier works instead of the 5f{k)/k ratios. Somewhat fortuitously, the numerical 
value of the perturbative difference in Eq.(19) with a Yukawa-type (A^ = Akf/Ti) screening 



5o^(M-5f(M = ^ 



(36) 



2kFj \ A2 

is not far, at least for Z = 1, from the precise value obtained from Eq.(26) with self- 
consistently determined phase shifts of an embedding problem. Qualitatively, the mistake 
one makes in choosing a linearly screened potential is "compensated" by the use of the 
first-order Born approximation; neither of these approximations are quantitatively accurate. 

More importantly, the careful numerical analysis of Ref. [46], performed by assuming 
higher (/ > 1) angular-momentum channels (i.e., using more ai parameters) to an estimation 
on TLS'-coupling, shows that (contrary to naive earlier expectations) there is not a series of 
Kondo-like effects corresponding to the increase of different orbital channels. This, numerics- 
based statement, is in harmony with the rule established at Eq. [24] in the present paper. 

A second-order Born approximation [valid for (Z/kp) < 1] for the scattering (transition) 
amplitude [f^"^^] with the above simple Yukawa potential gives 

2./'^'(,.M= V^W + ^,^^g^, (37) 

in atomic units, for simplicity. This approximation with a fixed, linear-response-based input 
potential would suggest an enhancement of the coupling. A more consistent (cons) treat- 
ment, in which the linear screening is also modified upto the second-order by a quadratic- 
response method, reduces^i this enhancement as follows 

2Z^A 

f^ici. k,) = f^^\q, k,) - xiTi^AM^T^^TW ^^'^ 

showing that fcons{q,kF) > V^^(g)/27r still holds for g 7^ 0, but in the forward {1!} = 0) 
limit fmns{q = 0, kp) = V^^{q = 0)/27r. This observation, which is based on selected (RPA) 
diagrams, heralds that care is needed when one uses a higher-order method in terms of a 
bare (in our case: a linearly screened) input potential in field-theoretic many-body attempts. 

A central local potential, for example Vps{r), gives rise to scattering of all orders of 
spherical harmonics even at the first-order Born level. On the other hand, a so-called 
separable potential^>^ for a given orbital (/) momentum causes scattering only for the 
given (/th) partial wave; in the case of all I, one speaks of a completely separable potential. 
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A correct determination of the corresponding channel-potentials, Vi{ki, /C2) = yiui{ki)ui{k2) , 
in this modeling could rest on experimental data or on a detailed, microscopic theory. 

It is important from the point of view of physical consistency, that the frequently applied 
tangent-method^ on multiple scattering effects gives^ a tan6o{kp) in terms of kp, u^ikp), 
and Vq of the s-channel potential. In the even simpler case with a contact Vo(r) ^ Vq5{y) 
auxiliary interaction (which gives a constant potential in momentum space) one has the 
remarkably simple form tan(5o = — TrpoVo for the s-wave phase shift. Note that by writing, 
formally, tan^o in the Ihs of Eq.(18), one can get this exact result in one-step from 



- 2mkF r , 2 
tanoo = — / drr 



6{r) 



4^,2--. foikFr). (39) 

The precise derivations of the exact result rest on much more involved calculations by us- 
ing real-space^ or momentum-space^ Schrodinger equation with Vo6{r). Application of a 
renormalized contact interaction, defined via — Vq ^ [V(^Po)] tan^o, in the perturbative rhs 
of Eq.(39) results in an identity. Clearly, this renormalized interaction could be used with 
unperturbed plane wave states to calculate the real reactance-matrix. The exact standing 
(st) wave solution [cf. Eq.(45)] of the tangent-method is the following 

(f)st{kF,r) = 1 . (40) 

cosoo kpr 

This is based, oncemore, on the principal-value Green's function. 

The standard logic to determine a value of Vq is based on the phase shift [5o(^)] of the 
real potential [Vo(r)], but the s-channel contact interaction^^ has an important limitation 
when we apply it to the screening problem of a charge Z. It was pointed out earlier— 
that this model cannot supply enough charge to shield the Coulomb field of the physically 
simplest impurity Z = 1. A formal requirement of 6o{kp) = 60 would result in Vq 00, 
when the self-consistent [5o(^f)] leading phase shift (for Z = 1) goes to 7r/2. At = 2.5 of 
the Wigner-Seitz density parameter, an about Vq ~ 20 value (in a.u.) is prescribed by this 
formal constraint. 

The auxiliary contact-potential was applied earlier within a local-phonon model of 
electron-phonon interactions^ in an opposite way, i.e., via a direct approximation [Vq ~ d) 
for it to characterize the (constant) momentum-space potential of an atomic displacement 
in a free electron gas. In this case tan^o could measure a renormalized effect of a prefixed 
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input Vo, beyond the conventional first-order Born approximation {Sq^^ = — vrpoVo)- Moti- 
vated by that work, the scahng equations of the two-level problem has been investigated^*^ 
in details in terms of an input 6q < 7r/2 variable in the so-called large-phase-shift case. The 
mentioned direct approximation was implemented^ recently as tan 60 = tt with a dimen- 
sionless matrix element {v^ ~ d) of the Coulomb-potential gradient of an embedded proton 
taken between the corresponding hydrogenic (Is and 2p) bound states, and weighted by an 
enhanced density of states due to an other model potential; cf. Ref. [55]. 

The scattering with contact interaction is isotrop, i.e., the scattering amplitude (and thus 
the diagonal transition T-matrix) does not depend on the scattering angle. The correspond- 
ing, dimensionless [kpFolkp)] exact scattering amplitude is 

kpFoikj.) = e^^» sin^o = - , . • (41) 

I — t tan do 

Its real part has the simple form of 

kp ReFoikp) = rrpo ^ ^ /° (42) 
1 + (ttpoVo)^ 

which could suggest, based on a proper reinterpret at ion of the rhs, an effective (contact) 
potential to calculations with unperturbed plane wave states. If the impurity scattering is 
purely local one can use a many-body Green-function-based method^ also to derive the 
exact result in Eq.(41), as the T-matrix depends only on energy and not on momentum; 
such a derivation shows the algebraic nature of the propagator method. 

For a single impurity embedded in an electron gas, and described by a regular Vo(r) 
potential, the electron Green's function can be written symbolically as 

G = Go + GoTGo, (43) 

using the standard T-matrix approach^. In our case the bare Green's function of a free 
electron of the ideal system with chemical potential // is GQ{Lj,'k) = [u; — (ek — /i) + iO]~^. 
Since the Fourier transform of a contact interaction is Vq, the equation for the corresponding 
T-matrix is particularly simple as it depends only on energy and not on momentum 

T{uj) = Vo + Vo J2Go{uJ,k)T{uj). (44) 
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One easily gets the T = Vo/(l — VoGq) solution, where the Gq^uj) quantity is the local 
propagator at the impurity site averaged in the momentum-vector and taken at the Fermi 
surface, |k| = kp and u = 0. Thus one has^ the simple Gq = —iirpQ relation, which results 
in Eq.(41) for the complex scattering amplitude; Trpo = mA;^/(27r/i^). Finally, the exact 
solution for the scattered wave, with Vo(r) = Vo6{r — Fq) and at A; = kp, is given by 



showing a limited applicability to a realistic enhancement at the impurity position. 

The main point in such a standard calculation is to show that the Green function method 
can, sometimes, demonstrate very general properties of many-body systems; here only the 
density of states of the unperturbed system comes into the final answer in Eq.(44). It is, on 
the other hand, not necessarily the easiest way of calculating beyond the contact-form for 
an impurity potential. The propagator, being essentially a function of two vector variables, 
is inevitably more complicated mathematically than a one-particle wave function. 

V. APPENDIX B: MULTIPLE SCATTERING OFF A LOCALIZED AV 

Until now the assumption has been used that the perturbation due to AV is weak, which 
is certainly the case for small d and valid in the parameter range has been discussed. The 
effect of this perturbation was considered within a distorted- wave Born approximation, i.e., 
via the first term of the distorted-wave Born series; the outgoing state of V is used for the 
final state of -|- AV. A precise treatment for arbitrary d, as we mentioned in Sec. Ill, 
should rest on a set of coupled radial equations of an axially symmetric potential. Such a 
treatment could model faithfully the d-dependence of coupling due to atom-displacements. 

For completeness, it is worthwile to discuss the multiple scattering in AV which is also 
makes connection to some related earlier works^i^i^. We outline, therefore, a Green- 
function-based analysis, supposing a convenient strongly localized potential-form and con- 
vergency of the corresponding distorted-wave Born series. Using a schematical notation, 
V + AV, the Dyson equation for the total Green's function G is as follows 




(45) 



G = Go +Go{V + AV) G 



(46) 
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where Go refers to the free electron. That equation can be solved in two steps. First, the 
summation is made only in V and an intermediate Green's function (corresponding to 
propagation in V) is defined symbolically as 

'''' - (*^' 

That can be expressed^^ by the exact independent-electron wave functions ip^lr), given 
by Eq.(6) in the presence of the self- consistent V in an electron gas as 

After performing the integration with respect to the direction of k, one gets 

G^(r, v\uj) = —^y [ dk e ^^'-^"^^^^'-^'^ . (49) 

where fkim = Ri{k,r)Yim{r)- In the usual (m = 0) approximation this equation is simplified 
(see Appendix A) and results in 

G^(r,r',t^ = 0) = -iiTpo ^ (pfc^,(r')v5fcj.«(r), (50) 

which is a generalization of the noninteracting Green function based on plane waves. 
Second, the term AV is taken into account, in analogy of Eq.(43), as 

^ = 1 - AVGV = + G^AK//G^ (51) 
where a nonlocal effective {ejf) term is introduced as 

In the general case Eq.(51) in real space is very complex, thus only the special case of a 
well-localized AV is treated; G^(0, 0, 0) = —inpoE^kp) where E^kp) is an enhancement due 
to V . Thus, the matrix elements of a strongly localized AV are calculated. 

The wave functions ^ki{'^) are regular (see. Fig. 2), thus they can be expanded in r and 
only the first corrections are kept {kpr <^ 1). With such constraints one gets 

d 

^kpi{^) = ^ifi ^kpi{^ = 0) + 5i^i r ■ —ipkpiir = 0). (53) 
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The matrix of AV can be calculated as e.g. 



= V>k,i=o{r = 0) 1 (^^^ ) r AV{r) dr. (54) 

Note that AVi^o = AVo^i. Furthermore, AVqo = as the volume integral of a dipolar AV{r) 
is zero, and AVii is 0{d'^k^). In such a way the matrix element of the effective potential 
can be [vrpo = kp/{2n)] determined as 

Therefore, for the physically realistic small d (where AV^ ~ d) the above-outlined treatment 
heralds a third-order change in d beyond the ci-linear term. Interestingly, the form in Eq.(55) 
is quite similar to Eq.(42). Finally, as we discussed in Appendix A, a convenient phase shift 
(do) could be introduced as Sq = — arctan(7rpoAVo,i). 
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